Controlled spin domain creation by phase separation 
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We demonstrate a method of controlled creation of spin domains in spin-1 antiferromagnetic Bose- 
Einstein condensates. The method exploits the phenomenon of phase separation of spin components 
in an external potential. By using an appropriate time dependent potential, a composition of spin 
domains can be created, as we demonstrate in the particular cases of a double well and a periodic 
potential. In contrast to other methods, which rely on spatially inhomogeneous magnetic fields, 
here the domain structure is completely determined by the optical fields, which makes the method 
versatile and reconfigurable. It allows for creation of domains of various sizes, with the spatial 
resolution limited by the spin healing length only. 

PACS numbers: 03.75.Mn, 03.75.Hh, 67.85.Bc, 67.85. Fg 



I. INTRODUCTION 

Bose-Einstein condensates with spin degrees of free- 
dom [l[ attracted in recent years great interest due 
to the unique possibility of exploring fundamental con- 
cepts of quantum mechanics in a remarkably controllable 
and tunable environment. The ability to generate spin 
squeezing and entanglement [2\ makes spinor Bose gases 
promising candidates for applications as quantum sim- 
ulators i3|, in quantum information jilj, and for precise 
measurements [5|. Moreover, spinor condensates were 
successfully used to recreate many of the phenomena of 
condensed matter physics in experiments displaying an 
unprecedented level of control over the quantum system. 
In particular, spin domains spin mixing Q, and 

spin vortices [10| were predicted and observed. 

The ability to create spin domains is a crucial com- 
ponent of applications including data storage and spin 
based logic implementation To date, domain struc- 
tures in Bose-Einstein condensates were created in a con- 
trollable way using inhomogeneous magnetic fields @ . In 
principle, these can be induced by magnetic coils, elec- 
tronic chips or permanent magnets [12]. However, both 
the spatial resolution of this method and the ability to 
produce arbitrarily shaped, reconfigurable spin domain 
structures is severely limited. 

In recent papers , [l^ , we considered the possibility 
of phase separation in the ground state of spinor con- 
densates. We demonstrated that this phenomenon can 
take place in both antiferromagnetic and ferromagnetic 
condensates in external potentials. In this paper, we ex- 
ploit the phenomenon of phase separation to propose a 
method of controlled creation of metastable spin domain 
states with a chosen spatial spin structure. It consists 
of applying an appropriately designed optical potential, 
which is subsequently slowly relaxed towards the desired 
final trapping potential. In contrast to the other method, 
relying on magnetic field gradients to separate spin com- 
ponents, our method uses more easily reconfigurable op- 
tical laser fields to shape the spatial structure. Impor- 
tantly, these structures are not "pinned" by local extrema 
of the magnetic field. An additional advantage is the high 



spatial resolution, which we demonstrate to be generally 
limited only by the spin healing length of the condensate. 

The paper is organized as follows. Section reviews 
the mean- field model of a spin-1 condensate in a homoge- 
neous magnetic field. In Section Hill we present numerical 
simulations of the spin domain creation process. In Sec- 
tion IIVI we discuss practical limitations of the proposed 
method using both analytical and numerical approach. 
Section IVl concludes the paper. 



II. MODEL 

We consider a spin-1 Bose-Einstein condensate in a ho- 
mogeneous magnetic field pointing along the z axis. We 
apply the mean field approximation, which describes di- 
lute Bose-Einstein condensates at zero temperature. In 
spinor condensates, ground states can substantially devi- 
ate from the mean field solutions even at zero tempera- 
ture due to large quantum fluctuations [^. However, it 
was shown that the introduction of a weak magnetic field 
restores the validity of the mean-field model. We start 
with the Hamiltonian H = Hg + Ha, 



H = 



E 



-,1 , Co 



(1) 

where t/j^ , ■(/^o, "0+ ^-re the wave functions of atoms in mag- 
netic sublevels mf = —1,0,-1-1, m is the atomic mass, 
V{v) is an external potential and ^ = ^77,^ =^ l^jP is 
the total atom density. The asymmetric (spin dependent) 
part of the Hamiltonian is given by 



C2 , 



(2) 



where Ej is the Zeeman energy shift for state tpj and the 
spin density is, 

F = {F.,,Fy,F,) = {^j^F,iP,^^Fy^,^P^F,^P) (3) 

where F^^y^z are the spin-1 matrices [l6l | and ■0 = 
(-0+, 7 0-)- The spin- independent and spin-dependent 



2 



interaction coefScients are given by cq = 47r?i^(2a2 + 
ao)/3m and C2 = AiTh?{a2 — ao)/3m, where as is the 
s-wave scattering length for colliding atoms with total 
spin S. The total number of atoms and the total magne- 
tization in the direction of the magnetic field 

N = j ndv, (4) 

j F,dr = J (n+-n_)dr, (5) 

are conserved quantities. The Zeeman energy shift for 
each of the components, Ej can be calculated using the 
Breit-Rabi formula Il7| 



(6) 



En — —-E] 



HPS 



where -Bhfs is the hyperfine energy splitting at zero mag- 
netic field, a = (gj — gi)^-BB j Eyiys, where /is is the 
Bohr magneton, gj and gj are the nuclear and electronic 
g-factors, and B is the magnetic field strength. The lin- 
ear part of the Zeeman effect gives rise to an overall shift 
of the energy, and so we can remove it with the transfor- 
mation 
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FIG. 1: (a) Ground state of a sodium condensate with fixed 
magnetization M — 0.4, number of atoms A'^ = 8.4 x 10'^, 
trap frequencies cjp = 27r x 1000 Hz and oj^ = 27r x 2.5 Hz, 
and the magnetic field strength B — 120 mG. The n+, no, 
and n components are depicted by dash-dotted, dashed, and 
dotted lines, respectively, and the solid lines show the total 
density, (b) The result of a slow introduction of a Gaussian 
shaped barrier with A = 3.3 x 10^^" J and w — 90/im during 
the time t = 5 s. 



in a ^^Rb condensate far from Feshbach resonances [T^. 
Consequently, changing the total density n requires much 
more energy than changing the relative populations of 
spin states Uj. In our considerations we will assume 
that the total atom density profile n{r,t) is close to the 
Thomas-Fermi profile for a given potential V{r,t). 



H ^ H + {N + M)E+/2+ {N - M)E-/2. (7) 

This transformation is equivalent to the removal of the 
Larmor precession of the spin vector around the z axis 
[isj . We thus consider only the effects of the quadratic 
Zeeman shift. For sufficiently weak magnetic field we 
can approximate it by SE = {E+ + E^ — 2Eq)/2 w 
a^£'HFs/16, which is positive for *^Rb and ^'^Na con- 
densates. 

The Hamiltonian ([1]) gives rise to the Gross-Pitaevskii 
equations describing the mean-field dynamics of the sys- 
tem 

dip 

'^^~gf ^[^ + C2{n± + no - n=p)] ^^^± + C2ipoip^ , (8) 
dip 

"^^-Qf = [C-5E + C2{n+ + n-)] ■00 + 2c2'(/'+'0_'0o ' 

where £ is given by £ = —h^V^ /2m + cqii + V{r). 

By comparing the kinetic energy with the inter- 
action energy, we can determine the healing length 
^ = 2'Kh/ \/2mcon and the spin healing length — 
2'Kh/ \/2niC2n. These quantities give the length scales 
of spatial variations in the condensate profile induced by 
the spin-independent or spin-dependent interactions, re- 
spectively. Analogously, we define the magnetic healing 
length as f s — 2T:h/ y/2mSE. 

In spin-1 condensates created to date, the ao and 02 
scattering lengths have similar magnitudes. The spin- 
dependent interaction coefficient C2 is then much smaller 
than its spin-independent counterpart Cq. For example, 
this ratio is about 1:30 in a ^^Na condensate and 1:220 



III. SPIN DOMAIN CREATION 

In Refs. [H, [l3| we demonstrated that in the pres- 
ence of a homogeneous magnetic field, a spin-imbalanced 
(magnetized) antiferromagnetic condensate is subject to 
phase separation in the ground state resulting in sepa- 
rate spin domains. The spin domains align according to 
a simple rule; the magnetized domains reside in the re- 
gions with the lowest density, while the unmagnetized 
domains fill the areas of high density. In this way, the 
nonlinear spin energy of the antiferromagnetic conden- 
sate is minimized for a given total magnetization. In 
this work, we will show that one can use this property to 
shape the spatial distribution of spin with an appropriate 
time varying external potential. 

We consider a quasi-ID sodium condensate trapped in 
an elongated harmonic potential described by the ID ver- 
sion of Eqs. ([5]) [ij] with rescaled interaction coefficients 
cJ^,C2^ — (m''^±)/(27r?i)co, C2, wherew_L is the transverse 
trapping frequency. The Fermi radius of the transverse 
trapping potential is smaller than the spin healing length, 
and the nonlinear energy scale is much smaller than the 
transverse trap energy scale, which allows us to reduce 
the problem to one spatial dimension [l^, [23|. To cre- 
ate a double well potential, we add a Gaussian potential 
barrier that can be realized eg. by an additional blue- 
detuned laser beam 

Viz, t) = imw^z^ -I- A{t)e~'-'/'^'"' . (9) 
In the absence of the barrier, and under certain exper- 
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FIG. 2: (a) Ground state of the double well potential with 
parameters as in Fig. [1] (b) Ground state after the removal 
of the Gaussian barrier, (c) Spin domain structure after the 
removal of the Gaussian barrier within t = 5s. (d)-(f) Same 
as (a)-(c), but for a weaker magnetic field B — 45 mG and 
magnetization M — 0.8. 



imental conditions [l3|, the ground state is character- 
ized by spatial separation of the m/ = and ruf — +1 
atoms (or my = —1, depending on the magnetization), 
see Fig.[lja). After switching the barrier on, the conden- 
sate splits in an effective longitudinal double well poten- 
tial into two symmetric parts. We note that the ground 
state now has the structure corresponding to Fig. ^a) . 
However, the result of a slow increase of the barrier on 
the timescale scale of several seconds results in a different 
state, where inner irif = +1 domains are absent on both 
sides of the barrier. The difference between Fig.[IJb) and 
Fig. [UJa) is a result of the strong spin immiscibility of 
the TO/ = 0, atoms that suppresses the tunnelling of 
rUf = +1 atoms through the mf — Q atoms. The final 
state is thus a metastable state, but practically stable on 
the relevant experimental timescales. 

A yet more intriguing effect can occur when the con- 
densate is prepared in the ground state of the double 
well potential (with the barrier on). Fig. [^Ja), and after 
gradually lowering the potential barrier. The condensate, 
instead of evolving towards the single-well ground state, 
Fig.IHb), creates another kind of a metastable state, with 
five alternating spin domains of my =0 and mf ~ +1 
atoms shown in Fig. ^c) . Again, we checked that these 
spin domains are perfectly stable on the timescales as 
long as several seconds due to suppressed tunneling. In 
principle, by applying several potential barriers with cer- 
tain parameters, it is possible to create a given number 
of stable spin domains of various sizes. 
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FIG. 3: (a) Ground state in a periodic potential with 
d = 140 ^im and A = 3.3 x 10"^° J. (b) The resuhing domain 
structure after switching the potential off during f = 5 s. 



We note that the above spin-domain creation scenario 
is possible only for domains which consist of immiscible 
atoms. The phase separation of to/ = and to/ = ±1 
atoms is a feature of antiferromagnetic spin-1 conden- 
sates in a relatively strong magnetic field [ij]. In the 
weaker field regime, the spin-imbalanced antiferromag- 
netic condensate generally consists of to/ = -1-1 and 
TO/ — —1 atoms, which are miscible. In result, the spin 
domains do not separate, and the slow process always 
results in a state close to the ground state of the system, 
as shown in Fig. [2Id)-(f). 

Furthermore, we consider a condensate placed in a pe- 
riodic optical potential 



V{z) ^ A{t) cos^{2Trz/d) , 



(10) 



with d being the lattice period. We simulate the dy- 
namics of the condensate with periodic boundary con- 
ditions, which can be experimentally realized eg. with a 
ring-shaped trapping geometry [2l[. Analogously as in 
the previous case, we start with a condensate with the 
optical lattice switched on. the density and spin pattern 
is shown in Fig.[31[a). By gradually decreasing the optical 
lattice strength, we arrive at a metastable state depicted 
in Fig. EKb), composed of multiple to/ = and to/ = -1-1 
spin domains. This state is also stable on the timescales 
of the order of seconds. We note that by modifying the 
lattice potential through introducing higher order Fourier 
components, we are also able to adjust the size of indi- 
vidual domains. 



IV. LIMITATIONS OF THE METHOD 

We now discuss some limitations of the above method 
of domain creation. As we mentioned above, these 
metastable states are stable on an experimentally rele- 
vant time scale due to the suppression of tunneling. How- 
ever, the tunneling rate depends on the system parame- 
ters such as the trap size and atom number. In Fig.[?]we 
show the result of adiabatic barrier removal for tighter 
harmonic traps with a;^ — 10 Hz (a) and ujz = 20 Hz 
(b). The atom number is here reduced proportionally to 
keep a similar density profile. Thus, an analogous domain 
structure as in Fig. ^c) is expected. Since the healing 
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length, determining the size of domain boundaries (see 
below) is almost the same in all cases due to a similar 
atom density profile, we expect the tunneling to be more 
significant in smaller systems. Indeed, the effect of the 
tunneling is clearly visible when comparing Fig.[2jc) with 
Fig.SKa) and (b). With the decreasing distance between 
the domains the state becomes unstable, see Fig. IH^a). 
In the case of the tightest trap. Fig. Hljb), the middle 
domain dissolves completely and the single-well ground 
state is obtained. 

The tunneling rate can be estimated analytically along 
the lines of Ref . Q , where the case of a spatially varying 
magnetic field potential was considered. Here, we apply a 
similar method to the case of a spatially varying external 
potential. First, we assume that the condensate separates 
into two components, one of them being ruf = and the 
other either mf = +1 or mf = —1, and neglect the 
influence of the third component. The potential energy 
for the ruf ~ i atoms is given by 



V^{z) = V{z) + giUi + giji 



ill) 



where j is the other present component. The interac- 
tion constants in an = 1 condensate are go = cq and 
g±i = go±i = Co + C2. We estimate the tunneling rate of 
ruf — atoms from the middle domain in a structure 
similar as the one shown in Fig. [2Ic) towards the edges 
of the condensate. The chemical potential of atoms in 
the middle domain is ^l^i = V{z) -I- gin+i(z) — const. 
We assume that the density of nif — component 
is negligible. Analogously, in the neighboring domain 
the nif — Q atoms have chemical potential equal to 
Mo — ^(-2^) + ffo'^oC-^)- The m/ — +\ atoms which 
tunnel through this domain feel the potential Vj^i{z) — 
V{z) -|- (701^0(2). Their density profile can be calculated 
in the WKB approximation 



dn^i{z) 



dz 



(12) 



Vmo - M+i + "-o(^)(.9oi - .90) 



where the decay of the density means that dn^i/dz > 
on the right side of the middle domain and dn^i/dz < 
on the left side. In the simplest approximation, uq^z) 
does not vary significantly and the decay is exponential. 



FIG. 5: Ground state of the double well potential for param- 
eters as in Fig. [If a) but with w — 9.4/im. 



Taking into account the pressure balance at the domain 
boundary, gono{zo)^ — gin+i{zQ)'^ , we can estimate the 
decay as n ~ exp(— az) with a = 2A/mc2no(zo)/ft — 
2(3/2)7^^-1^ where is the spin healing length. We can 
now write down the formula for the tunneling rate in the 
metastable state 



dN 
~dt 



(13) 



where Az is the distance between the domains and 7 is 
the attempt rate 7 = 2(~'^/^^no(2;o)(^'s), with (vg) being 
the speed of sound averaged over the transverse density 
profile [7,]. This indicates that the tunneling is greatly 
suppressed when the distance between the domains is 
larger than the spin healing length, which determines the 
spin domain boundary width through the coefficient a. 
In Table [J we show the ratio of the calculated amount of 
atoms that tunnel through to the total number of atoms 
in the middle domain. In the case of the tightest trap, the 
tunneling over the time of the evolution would exceed the 
number of available atoms, which explains the absence of 
a metastable state. 



Trap frequency 


-/Vtuniiclod /A^total 


2.5 Hz 


10-^2 


10 Hz 


0.05 


20 Hz 


2 



TABLE I: The ratio of the number of atoms that can tunnel 
through the barrier according to Eq. (|13|) within t = 5 s to the 
number of available atoms in the middle m/ = +1 domain. 

The second limitation of the method is related to the 
minimal size of a domain. For illustration, in Fig. [5] we 
present the condensate ground state in the double well 
potential as in Fig. [2{a), but created with a very narrow 
Gaussian barrier. In this case the inner domains of mj — 
+1 atoms are absent, and consequently the slow decrease 
of the barrier height leads to the ground state. Fig. [2Kb), 
instead of the metastable state. Fig. HJc). The reason 
for the absence of the domains is the very steep slope 
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of the density profile close to the barrier. According to 
Ref . (14] , in the local density approximation (LDA) , the 
spin state at a given point is determined by the value of 
the atom density at this point. The mj — +1 atoms can 
reside only in a very narrow region on the density slope 
where it is smaller than a certain value. If this region is 
smaller than twice the spin healing length, the domain 
cannot form. The approximate size of this region can 
be estimated as twice the distance between the point of 
maximum density n,nax and half the maximum density, 



2^/21n2 
Az « — =w . 

In^— 



(14) 



We confirmed through systematic numerical simulations 
with varying barrier size that with a good accuracy the 
inner domains appear when Az becomes larger than 2^s. 



netic condensate in a homogeneous magnetic field. The 
method exploits the phase separation of spin components 
in an external potential. In contrast to other methods, 
which rely on the spatially varying magnetic fields, the 
domain structure is here completely determined by op- 
tical fields, which makes this method more versatile and 
reconfigurable. Additionally, the method allows for cre- 
ation of domains of various sizes, with spatial resolution 
limited by the spin healing length only. 
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